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Abstract 

We consider simple hydrodynamical models of galactic dark matter in which 
the galactic halo is a self-gravitating and self-interacting gas that dominates the 
dynamics of the galaxy. Modeling this halo as a sphericaly symmetric and static 
perfect fluid satisfying the field equations of General Relativity, visible barionic mat¬ 
ter can be treated as “test particles” in the geometry of this field. We show that 
the assumption of an empirical “universal rotation curve” that fits a wide variety 
of galaxies is compatible, under suitable approximations, with state variables char¬ 
acteristic of a non-relativistic Maxwell-Boltzmann gas that becomes an isothermal 
sphere in the Newtonian limit. Consistency criteria lead to a minimal bound for 
particle masses in the range 30 eV < m < 60 eV and to a constraint between the 
central temperature and the particles mass. The allowed mass range includes popu¬ 
lar supersymmetric particle candidates, such as the neutralino, axino and gravitino, 
as well as lighter particles (m ~ keV) proposed by numerical N-body simulations 
associated with self-interactive CDM and WDM structure formation theories. 
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1 Introduction 


The presence of large amounts of dark matter at the galactic lengthscale is already an 
established fact. Assuming that this dark matter is a gas (or gas mixture) of various 
particles species, the established classification criteria labels possible dark matter forms 
as “cold” or “hot”, depending on the relativistic or non-relativistic nature of the particles’ 
energetic spectrum at their decoupling from the cosmic mixture |^, |0, [Q. Hot dark mat¬ 
ter (HDM) scenarios seem to be incompatible with current theories of structure formation 
and thus, are not favoured dark matter candidates |^, [Q, |Q. Cold dark matter (CDM), 
usualy examined within a newtonian framework, can be considered as non-interactive (a 
self gravitating gas of collisionless particles) or self-interactive [Q . CDM models are often 
developed in terms of n-body numerical simmulations 0, [0, [§, 0. Non-interactive 
CDM models present the following discrepancies with observations at the galactic scale 
0, O: (a) the “substructure problem” related to excess clustering on sub-galactic 
scales, (b) the “cusp problem” characterized by a monotonic increase of density towards 
the center of halos, leading to excessively concentrated cores. These problems appear in 
the more recent numerical simulations (see 0], [§, 0). In order to deal with these prob¬ 
lems, the possibility of self-interacting dark matter has been considered, so that nonzero 
pressure or thermal effects can emerge, thus leading to self-interactive models of CDM 
(z.e. SCDM) lil, [|T|, 


14| , ||15|| , fTql and “warm” dark matter (WDM) models 


that challenges the duality CDM vs. HDM. Other proposed dark matter sources consist 
replacing the gas of particles approach by scalar fields |0|, and even more “exotic” 
sources 0]- 

Whether based on SCDM or WDM, current theories of structure formation point 
towards dark matter characterized by particles having a mass of the order of at least keV’s 
(see 


T^-|2^), thus suggesting that massive but light particles, such as electron neutrinos 


and axions (see Table 1), should be eliminated as primary dark matter candidates (though 
there is no reason to assume that these particles would be absent in galactic halos). Of all 
possible particle candidates (denoted as WIMP’s: weakly interactive massive particles) 
complying with the required mass value of relique gases, only the massive neutrinos (the 
muon or tan neutrinos), have been detected, whereas other WIMPS (gravitino, sterile 

121 , 


neutrino, axino, etc.) are speculative. See ^6 


81, IM and Table 1 for a list of 


candidate particles and appropriate references. 

Even if the dynamics of visible matter in galaxies can be described succesfuly with 
Newtonian gravity, we believe that General Relativity is an appropriate framework for 
understanding basic features of galactic dark matter, a gravitational field source whose 
precise physical nature still remains an open question. If the results obtained with GR co¬ 
incide with Newtonian results, then there is no harm done from a pragmatic calculations- 
oriented point of view. However, from a formal-theoretical approach, we believe it is 
beneficial to broaden the scope of the study of galactic dynamics by incorporating it 
to a more general gravitational theory. In particular, in this paper we aim at testing 
the compatibility between observed galactic rotation curves and simple thermodynamical 
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assumptions under the framework of GR. 

Since dark matter halo probably constitutes the overwhelming mayority (90 %) of the 
galactic mass, an alternative approach to numerical simulations and newtonian hydro¬ 
dynamics follows by a general relativistic model describing the gravitational held of the 
galaxy as a spacetime geometry generated by the the dark matter halo (as a self gravi¬ 
tating gas), hence visible matter becomes test particles that evolve along stable geodesics 
of this spacetime. There is strong empiric evidence that the radial prohle of rotational 
velocities (“rotation curves”) in most galaxies roughly hts a “universal rotation curve” 
(URC) [^, IQ. This URC is characterized by a “battening” effect whereby rotation ve¬ 
locities tend to a constant “terminal” velocity whose value depends on the type of galaxy 
(between 125 km/sec and 250 km/sec). The profile of rotation curves identihes two main 
contributions of galactic matter: visible matter (the disk), showing a keplerian decay, and 
dark matter (the halo), explaining the battening ebect. This kinematic evidence might 
allow us to determine (at least partialy) the geometry of the spacetime associated with a 
self gravitating galaxy. In other words, our approach somehow inverts the standard initial 
value procedure in general relativistic hydrodynamics: instead of prescribing initial data 
based on physicaly motivated sources and then bnd the geometry of spacetime and the 
trayectories of test particles after solving Einstein’s equations, we provide brst constraints 
on the geometry of spacetime (from symmetry criteria and empirical kinematic data) and 
then bnd, with the help of the held equations, the corresponding momentum-energy ten¬ 
sor of the sources. This approach to galactic dark matter has been used in connection to 
scalar belds fQ . 


Bearing in mind that the dark matter halo overwhelmingly dominates the galactic 
matter content (at least in the halo region), we shall assume that the galactic halo (as a 
self gravitating gas) is the unique matter source of the galactic spacetime. Visible mat¬ 
ter becomes then test observers that follow stable circular geodesic orbits (the galactic 
rotation curves) of this spacetime. Following the “inverse” approach described above, 
we propose to use the empiric law governing the form of the URC for the galactic halo 
p0| and IQ) in order to make specibc asertions on the nature of the sources of the 


see 


galactic spacetime. Considering the self-gravitating galactic halo gas to be self-interactive 
(instead of colissionless matter or a scalar held), we aim at verifying if the assumption 
of the URC proble for the rotation velocity of geodesic observers (rotation curves) is 
compatible with the assumption that the galactic halo gas is a simple self-gravitating 
and self-interactive gas in thermodynamical equilibrium. For this purpose, we consider 
the galactic halo to be a spacetime characterized as: (a) sphericaly symmetric, (b) its 
energy-momentum tensor is that of a perfect buid satisfying the equation of state of an 

m 


equilibrium Maxwell-Boltzmann gas in its non-relativistic limit 


Assumption (a) 


is supported by observations in the halo of galaxies, while (b) is the central hypothesis in 
the present work. Assumption (b) requires spacetime to be stationary (static if rotation 
vanishes) and leads to a law relating temperature gradients with the 4-acceleration (Tol- 
man Law). Then, since we are assuming the validity of the empiric URC for the galactic 
halo, we need to cast the held equations and the conditions imposed by the thermodynam- 
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ics in terms of this rotation velocity [i.e. the velocity of test particles in stable geodesic 
orbits), considered now as a dynamical variable. Using the URC empiric law as an ansatz 
for this velocity immediately leads to expressions for the state variables that are (un¬ 
der suitable approximations) consistent with the thermodynamics of the non-relativistic 
Maxwell-Boltzmann ideal gas. From these expressions and bearing in mind numerical 
estimates of the empiric parameters appearing in the URC (“terminal” rotation velocities 
and the “core radius”), we obtain: (1) a constraint on the ratio of the particles mass 
to temperature for this gas, (2) the criterion of applicability of the Maxwell-Boltzmann 
distribution {i.e. the non-degeneracy criterion) leading to a minimal bound of about 
30 to 60 eV for the mass of the gas particles. Therefore, the assumption of a Maxwell- 
Boltzmann gas (SCDM or WDM) model for the galactic halo leads to an acceptable value 
for the particle’s mass lying in the range m > 0.5 keV. We provide in Table 1 a list of par¬ 
ticle candidates that could be accomodated according to the criteria (1) and (2) above, 
namely: neutralino, photino, light gravitino, sterile neutrino, dilaton, axino, majoron, 
mirror neutrino and possibly standard massive neutrinos. As mentioned previously, this 
mass range is compatible with predictions of current work based on SCDM and WDM 
structure formation models. We hnd it interesting to remark that barions and electrons 
comply with the criterion (2) above, but (1) would imply gas temperatures of the order 
of 10^ — 10® K. A gas of barions or electrons at such temperatures would certainly not 
be “dark”. HDM or WDM models based on less massive particles, like the electron neu¬ 
trino, remain outside the scope of the present work, since these particles might require 
assuming a fuly relativistic Maxwell-Boltzmann gas or a degenerate gas (possibly rela¬ 
tivistic) complying with a Fermi-Dirac or Bose-Einstein statistics. The axion, as well as 
other non-thermal relique sources, are also outside the scope of this paper and their study 
requires a different approach. 

The paper is composed as follows. In the next section we present the held equations 
for a static sphericaly symmetric spacetime with a perfect huid source. We provide in 
section 3 a review of the thermodynamics of an equilibrium Maxwell-Boltzmann gas. 
Then, in section 4, we re-write the held equations in terms of the orbital velocity of stable 
circular geodesics and then assume for this velocity the empiric ansatz given by the URC. 
This leads to forms of the state variables that will be compatible, under suitable series 
approximations, with the thermodynamics of the Maxwell-Boltzmann gas. Putting all 
these results together, we discuss in the last section the possible ranges for the mass of 
the particles of the dark halo gas and suggest future lines of research. 

2 Field Equations 

As mentioned above, we consider the halo to be spherically symmetric and is the de¬ 
termining component of the stress energy tensor determining the geometry. Thus, we 
consider the line element of an spherically symmetric space time: 
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+ (^de^ + sin^ e dif^) . 


( 2 . 1 ) 


ds‘^ 


-B^{r)c^dt^ + 


dr'^ 

1 — 2 M(r)/r 


Assuming as the source of such line element a static perfect fluid momentum energy tensor 
= {p + p) with 6°' ct, we obtain the following held equations 


GU = kp = 


2M' 




2 / 2M\ B' 2M 

G r = t^P = - 1 -- 
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where k = SirG/c'^ and a prime denotes derivative with respect to r. As mentioned in the 
Introduction (see also and [^), it is possible (by working with Einstein’s equations 
backwards) to impose constrains on the geometry of the spacetime that yield valuable 
information regarding the type of matter sources curving such spacetime. We apply this 
reasoning to the observed velocity prohle of stars orbiting around a galaxy, considered 
as test particles moving in stable circular geodesics. Knowing the specihc form of the 
velocity prohles around these geodesics should allow us to infere at least basic features 
on the nature sources producing the galactic held. For a sphericaly symmetric spacetime 
( 0 ). the energy, E, and the angular momentum, L, are conserved quantities for any 
particle moving in a geodesic, hence the geodesic equation for the radial motion has the 
following form: 


1 L2 \ 

B 2 (^^2 ^2 j 


0 , 


(2.5) 


where a dot denotes derivative with respect to the affine parameter of the geodesic. For a 
test particle to be in stable circular geodesic motion in any static spherically symmetric 
space time, its energy and angular momentum must satisfy: 


E^ 


53 

(2,6) 

B-rB' ’ 

r3 

(2.7) 

B-rB' ■ 


The tangential velocity of these test particles, V (r), can also be expressed in terms of the 
metric coefficients as: 
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( 2 . 8 ) 


o2 


= v^(r) = 


rB' 


B 

In previous work (see [^, [^] and |^) this tangential velocity was assumed to be 
constant along the full domain of the solution. In the present paper we consider v = v{r) 
and eliminate the metric coefficient, B{r), and its derivatives in (|2.2|) - (p.4|) in terms of the 
tangential velocity, n^(r), as given by Eq. (|2.8|) . This leads to a form for the held equations 
in which v{r) becomes a dynamical variable replacing B{r). 

After some algebraic manipulation, equation (^) becomes the following constraint 
relating the metric function M and the tangential velocity: 


,(—3 — 5v^ + Avv'r + 2v 
M' + ^^-r-- 


r (1 + 


1 + 


= u. 


(2.9) 


Substitution of this last equation into equations (|2.2|) and (p.3|) provides the following 
expressions for the density and the pressure of the huid in terms of M, v and v'\ 


up 
K p 


2 —M — 2 Mv"^ + v‘^r 


[—8 vv' r — 2 (2 + 1) — 3)] M + 4 r‘^vv' + 2v‘^ (—2 + r 

(1 + 


( 2 . 10 ) 


( 2 . 11 ) 


It is remarkable to see how the replacement of i? by n considerably simplihes the held 
equations for a general static and sphericaly symmetric held with a perfect huid source. 
Writing the held equations in terms of the orbital velocity, v, provides a useful insight 
into how an (in principle) observable quantity relates to spacetime curvature and with 
physical quantities (state variables) which characterize the source of spacetime. Thus, 
given an empirical functional form for v{r) (a rotation “prohle” for test particles), we can 
obtain M(r) by integrating the constraint (|2.9|) and thus, we arrive to fully determined 
forms of p and p in (|2.11|) . 


3 Thermodynamics 


We aim at verifying if the empiric laws associated with galactic rotation curves (hence, as¬ 
sociated with v{r)) can be compatible with matter sources that satisfy basic physical con¬ 
siderations and principles. Since galactic dark matter is, most probably, non-relativistic 
and the assumption of a perfect huid source for the static metric ( p.l|) points to an equi¬ 
librium conhguration, it is tempting to verify if p and p associated with a given empirical 
form for v{r) correspond to state variables characteristic of simple, non-relativistic sys¬ 
tems in thermodynamically equilibrium, such as a suitable ideal gas in its non-relativistic 
limit m, 0. 
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If we assume that the self gravitating ideal “dark” gas exists in physical conditions 
far from those in which the quantum properties of the gas particles are relevant, we 
would be demanding that these particles comply with Maxwell-Boltzmann (MB) statistics. 
Following 1^ 


the “non-degeneracy” condition that justihes an MB distribution is given 


by 


mk^ Tf/'^ 

where n, T, fi and kg are, respectively, the particle number density, absolute temperature, 
Planck’s and Boltzmann’s constants. If the constraint (|3.1|) holds and we further assume 
thermodynamical equilibrium and non-relativistic conditions, the ideal dark gas must 
satisfy the equation of state of a non-relativistic monatomic ideal gas 


3 

p = mc^ n + -nkgT, p = nkgT, (3.2) 

whose macroscopic state variables can be obtained from a MB distribution function under 
an equilibrium Kinetic theory approach (the non-relativistic and non-degenerate limit of 
the Jiittner distribution) ||^. An equilibrium MB distribution restricts the geometry of 
spacetime [Q, resulting in the existence of a timelike Killing vector field (5°" = (5 u°‘, where 
(3 = mc^/kgT, as well as the following relation (Tolman’s law) between the 4-acceleration 
and the temperature gradient 


Un 


+ h^(lnT),, = 0, = + 


leading to 


B' 

'b 


^ = 0 
T 


oc 


B 


-1 


(3.3) 


(3.4) 


The particle number density n trivialy satisfies the conservation law ;a = 0 where 
= nn“, thus the number of dark particles is conserved. Notice that given (|2.11|) , 
the equation of state (|3.2|) and the temperature from the Tolman law ( |3.4| ), we have two 
different expressions for n 


n = 


kgT 


n = 


mc^ 


oc pB, 


P-^P 


The quantity mc^ n in (p.6|) follows directly from equations (p.l 


(3.5) 

(3.6) 


Kmc n = 


[—8 vv' r -f- (n^ -|- 9) (2 -|- 1)] M -|- 4 vv' {7 + v^) r 

(1 -I- v^) 


(3.7) 


while n in ( |3.5| ) also follows from p in ( |2.11| ) with B oc exp[/ {y‘^/r)dr\. Consistency 
requires that (|3.5|) and ( p.6|) yield the same expression for n. 
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4 Dark fluid hydrodynamics 


So far we have expressed the held equations and the thermodynamics of the huid source 
in terms of the tangential velocities v{r) and the effective mass-energy M{r). If the dark 
matter component dominates the dynamics of the huid, we can ignore the contribution 
from visible matter (barions) and assume that the matter source is made exclusively 
of this dark matter component. A useful strategy to follow is then to prescribe, as a 
functional form for u(r), the empiric ansatzoi the radial prohle of tangential velocities of 
dark matter halo obtained from the “universal rotation curve” (URC) that roughly hts 
observed galactic rotation curves |^, This empirical form can be used in equations 
( ^■91) and ( |2.11|) . The function M(r) follows by solving the constraint (|2.9|) , subjected to 
the appropriate boundary conditions. The obtained M together with the prescribed v 
yield fully determined forms for p, p and the metric coefficient B{r). Next we can verify 
compatibility of n obtained from (|3.5| ) and (|3.6|). Finaly, we should be able to estimate 
the ratio (3 = mc^/k^T which in turn, from estimations of masses from particle physics, 
leads to an estimation of the temperature of the dark gas in terms of the particles’ mass. 

We shall assume for the empiric dark halo rotation velocity law given by Persic and 
Salucci [BDi 


m 


= 


Vq 

a? -\- 


X = 


opt 


(4.1) 


where r^pt is the “optical radius” containing 83 % of the galactic luminosity, whereas the 
empiric parameters a and uq, respectively, the ratio of “halo core radius” to and the 
“terminal” rotation velocity, depend on the galactic luminosity. For spiral galaxies we 
have: Vq = — /d*)(l + a3), where u^pt = n(ropt) and the best £t to rotation curves is 

obtained for: a = 1.5 and /3* = 0.72 -|- 0.441og]^o(-^/-^*)) where L* = lO^^'^L©. 

The range of these parameters for spiral galaxies is 125 km/sec < uo < 250 km/sec and 
0.6 < a < 2.3 ^ The field equation (p.9|) becomes the following linear first order ODE 


dM (1-1-2 vo'^) {vo‘^ — 3) {vg'^ + Q)x‘^ — 3 

dx X (6^ -1- x^) (a^ + x'^ + vo‘^x‘^) 

(2-^0^) Topt ^ 

(a^ -I- (a^ + x"^ + Vo‘^x‘^) 


the state variables (p.ll|) become 


(4.2) 


1 _ [(1 -F 2 vq'^) (3 - Vq'^) x"^ + {vp'^ + Q) x"^ + 3a^] M - vq"^ (2 - vp'^) x^ r„pt 

2 ^ {a'^ + x"^) [a? -I- (1 + vp'^) x"^] x^ 


1 _ -[a^ + {1 + 2 vp"^) x"^] M -f vp'^x^ r„pt 

2^^ (a^-t-a;2) rfpj 

^The URC given by CD fits also elliptic and irregular galaxies. 


(4.4) 
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while the metric coefficient B takes the form 


B 


1 + 



T = Tr 


1 + 


-vll2 


(4.5) 


SO that T is the temperature complying with the Tolman law and T^. = T(0). Given a 
solution M = M{x) of (|4.2|) , all state variables become determined as functions of x and 
Vq. The solution of ( |4.2| ) is the following quadrature 


(ug — 2 )(a^ + ^ovqX^ r + {1 + Vq)x‘^]^^ xdx 

[a2 + (1 + ^ (a^ + x^f-^’o 


(4.6) 


where we have set an integration constant to zero in order to comply with the consistency 
requirement that Uq = 0 implies flat spacetime (i? = 1, M = 0). Since the velocities of 
rotation curves are newtonian, Vq c (typical values are Vq/c ~ 0.5 x 10“^), instead of 
evaluating (|4.6|) we will expand this quadrature around vq/c (in order to keep the notation 
simple, we write Vq instead of Vq/c). This yields 


M = 


X3 y. 


opt 2 

vf 


2 I 2 


bx^ + 2a^ 2 + 3a^ 4 g 


2(a2 + a;2; 


■ + 


2(a2 + a;2)2 




(4.7) 


we obtain the expanded forms of p and p by inserting ( [4.6|) into (^4.3| ) and ( [4.4| ) and then 
expanding around uq, leading to 




2(3a2 + x^) 




(a2 + a;2)2 « 

2a2 + a;2 


2 5x + 23a X + 6a 4 . g. 

--Ug" + OivI), 


(a2 + a;2)3 


4 2a:^ + 7a2a;2 + 3a^ g g 

Vn --^77^-Ug + O(no), 


(a2 + x2)2 0 (a2 + a;2)3 

while the expanded form for T follows from ( [4.5|) 


T = Tr 


1 - i In (1 + vl + i lii^ (1 + x^) n)! - 0 (v‘) 


(4.8) 

(4.9) 

(4.10) 


In order to compare n obtained from ( p.5|) and (|3.6|) , we substitute (|4.1|) and (|4.6|) into 
( p.7|) and expand around vq, leading to 


nr , = 

opt 


nmc^ 


2(3a2 + 3)2) 
(a2 + x^) 


+ 55a2a;2 + 13a^ 4 


2^ ^0 


2(a2 + 0)2)3 


ti„ + o(!>; 


(4.11) 


while n in (|3.5|) follows by substituting ( |4.6|) into (|4.9|) , using T from ( |4.5| ) and then 
expending around Vq. This yields 
































nr . = — ; —— 
°Pt tc h T 

tv tv ^ -L Q 


2a^ + 


(a^ + x'^) 


2\2 ^0 




(2a^ + x‘^){a^ + x^) In (1 + x^) — 2{a? + 2x^)(3a^ + x^) 

2(a^ + x^)^ 


t.» + 0(1.^ 


(4.12) 


Since xq/c -C 1, a reasonable approximation is obtained if the leading terms of n from 
( [4.11| ) and ( [4.12|) coincide. By looking at these equations, it is evident that this consistency 
requirement implies 


(4.13) 


1 2 3 

-mvn = -k„ic 

2 ° 2 

where xq denotes a velocity (km/sec) and not the adimensional ratio xq/c. Equation (|4.13|) 
is analogous to the condition that dehnes the so-called “virialized temperature” in the 
context of cooling of a baryon gas, though in the present approach such a temperature 
corresponds to the dark matter gas (see section 17.3 of 0). Since higher order terms 
in Xq/c have a minor contribution, the two forms of n are approximately equal. This 
is shown in Figure 1 displaying the adimensional quantity nmc^nrl^^ from ( 4.11 ) and 
( [4.12| ) as functions of x for typical values Xq = 200 km/sec, a = 1 and eliminating 
with (|4.13|) . Equation ( [4.6|) shows how “flattened” rotation curves, as obtained from the 
empiric form (|4.1|) , lead to M oc for r 0 and M oc r for large r. Equations ( |4.3D to 
( [4.13| ) represent a relativistic generalization of the “isothermal sphere” that follows as the 
newtonian limit of an ideal Maxwell-Boltzmann characterized by p ~ mc^n,p <C p and 
T Tf.. In fact, using newtonian hydrodynamics we would have obtained only the leading 
terms of equations (|4.3|) to (|4.13|) . It is still interesting to hnd out that the isothermal 
sphere can be obtained from General Relativity in the limit xq/c <C 1 by demanding 
that rotation curves have a form like (H). The total mass of the galactic halo, usualy 
given as M evaluated at the radius r = r 2 oo (the radius at which p is 200 times the mean 
cosmic density). Assuming this density to be ~ 10“^® gm/cm^ together with typical values 
xo = 200 km/sec and a = 1 yields r 2 oo ~ 150 kpc. Evaluating M at this values yields 
about 10^^ Mq, while M evaluated at a typical Vopt = 15 kpc leads to about 10^^ Mq, an 
order of magnitude larger than the galactic mass due to visible matter. These values are 
consistent with Refs. and . 


5 Discussion. 

So far we have found a reasonable approximation for galactic dark matter to be described 
by a self gravitating Maxwell-Boltzmann gas, under the assumption of the empiric rotation 
velocity law €3). The following consistency relations emerge from equations 
and (|4.13|) 


(0T),(02) 
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3v^ 


AttG m ' 


Tr 


mv. 
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3 fc. 


(5.1) 


hence, bearing in mind that n < ric and T si T^, the condition (|3.1| ) for the validity of the 
MB distribntion can be written as 


n h'^ 


< 


Ur 


(mfc^T)3/2 - 

Inserting (|5.1|) into (|5.2| ) yields the condition 


< 1 , 


(5.2) 


m ^ 


3^/2 


47rG a2 no 


1/4 


(5.3) 


a criteria of aplicability of the MB distribution (non-degeneracy) that is entirely given in 
terms of m, the fundamental constants G, h and the empiric parameters Vq and b = 

(the “terminal” rotation velocity and the “core radius”). For dark matter dominated 
galaxies (spiral and low surface brightness (LSB)) these parameters have a small 
variation range: Ri ISkpc , 1 kpc <b< 5 kpc and 125 km/sec < Uq < 250 km/sec, 

the constraint (p.3D does provide a tight estimate of the minimal value for the mass of the 
particles under the assumption that these particles form a self gravitating ideal dark gas 
complying with MB statistics. As shown in Figure 2, this minimal value lies between 30 
and 60 eV, thus implying that appropriate particle candidates must have a much larger 
mass than this range of values. This minimal bound excludes, for instance, light mass 
particles such as the electron neutrino (rrij,^ < 2.2 eV) or the axion (m^ ~ 10“® eV). The 
currently accepted estimations of cosmological bounds on the sum of masses for the three 
active neutrino species is about 24 eV [^, a value that would apparently rule out all 
neutrino flavours. However, recent estimations of these cosmological bounds have raised 
this sum to about 1 keV [Q, hence more massive neutrinos could also be accomodated 
as dark matter particle candidates. Estimates of masses of various particle candidates are 
displayed in Table 1. 

Since T T^, the consistency condition (|4.13| ) provides the following constraint on 
the temperature and particles mass of the dark gas 


m 

C 



0.4 X 10^ 


T ’ 


(5.4) 


where we have taken uq = 200 km/sec ^ Considering in (|5.4|) the minimal mass range 
that follows from ( ^.3|) , we would obtain gas temperatures consistent with the assumed 
typical temperatures of relic gases: Tc ~ 2 to 4 K. However, as long as we do not have 
more information on the interaction and physical properties of various particle candidates, 
we cannot rule out a given large mass value on the grounds that the corresponding gas 

®The variation of vq in the observed ranges for spiral galaxies does not alter significanly the numerical 
value in the rhs of Eq. (^). 
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temperature could be too high. However, if we assume that the ideal dark gas is made of 
electrons or barions, so that m = rrip or m = rrie, then condition (|5.3|) for applicability of 
the MB distribution is certainly satished, but (|5.4|) implies a temperature of the order of 
Tc ~ 10^ K for electrons and T^. 10® K for barions ! Obviuosly, barions or electrons at 

such a high temperatures would certainly not remain unobservably “dark”. We can rule 
them out, but we cannot rule out more massive particles (in the range of 1-100 GeV’s) 
characterized by weak interaction even if the gas temperature is in the range of 10® — 10® 
K. Figure 3 illustrates, for various particle candidates, the relation between Tc and m 
contained in (^). The main novelty of the present paper is the fact that it is based on a 
general relativistic hydrodynamics, as opposed to numerical simulations 0-0) newtonian 
or Kinetic Theory perturbative approaches (see ||I^-||22||). 


Finaly, the fact that under the assumption of MB distribution, we have obtained 
a minimal mass on the range 30 eV < m < 60 eV that seems to discriminate against 
thermal relique gases composed by lighter particles (electron neutrino, etc) coincides with 
the fact that these particle candidates tend to be ruled out because of their inability 
to produce sufficient matter clustering 0, 0. In spite of these arguments, if either of 
these particles constitute a self gravitating gases accounting for a galactic halo it would 
be inconsistent to model such a gas as SCDM in the context of a classical ideal gas 
complying with MB statistics. It would be necessary to examine these cases as either 
HDM or WDM, by using either a relativistic MB distribution (very light particles can 
be relativistic even at low temperatures) and/or a distribution that takes into account 
(depending on the particle) Fermi-Dirac or Bose-Einstein statistics. Non-thermal axions 
are very light particles (m ~ 10“® eV), however this type of relique source cannot be 
treated as a Maxwell-Boltzmann gas, thus the lower mass limit that we have obtained 
does not apply. This and other non-thermal sources [^, require a wholy different 
approach. These studies will be undertaken in future papers. 
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SCDM/WDM 
Light Candidates 

mass in keV 

References 

Light Gravitino 

~ 0.5 

[0 


-0.75-1.5 

H, 113 


-2.6-5 

38 

Sterile Neutrino 

< 40 

s 


1 - 100 

10. 

Standard Neutrinos 

~ 1 

41, ia 

Dilaton 

- 0.5 

li 

Light Axino 

- 100 

H] 

Majoron 

~ 1 

Ml M 

Mirror Neutrinos 

~ 1 


CDM 

Heavy Candidates 

mass in GeV 

References 

Neutralino 

> 32.3 

s 


> 46 

1, 

Axino 

- 10 

52, 53 

Gravitino 

- 100 

|4 


Table 1: Particle candidates for a MB Dark Matter gas. 


12 








































































































Kmc^r 2 n 

op 



Figure 1: Comparison of n obtained from (^4.11|) and ( [4.12|) . This plot displays the adi- 
mensional quantity Kmc^r'^p n, as a function of x, obtained from truncating the right hand 
sides of ( |4.11| ) (solid curve) and ( [4.121 ) (dotted curve with crosses) up to fourth order in 
Vq/c and assuming the consistency condition ( |4-13|) . Since we are using an empiric law 
for observed galactic rotation curves (the URC given by ([4.1|)), the fact that these two 
expressions for n are so close to each other provides an empirical justihcation for the 
compatibility between MB distribution and these rotation curves. 


13 




















Figure 2: Minimal mass for which the Maxwell-Boltzmann distribution is applicable. 
This graph displays m (in keV’s) as a function of Vq/c and b = ax, respectively, the ter¬ 
minal velocity and ‘halo core radius’ associated with the URC given in (|4.1|) . Assuming 
typical ranges for spiral galaxies: 125 km/sec < vo < 250 km/sec and Ikpc <b< 5kpc, 
we obtain masses in the range of 30 eV < m < 60 eV that follow from the right hand side 
of the relation (^.31) , providing the criterion for applicability of the Maxwell-Bolzmann 
distribution. Dark matter particle candidates complying with an MB distribution must 
have much larger mass than the plotted values 30 eV < m < 60 eV. 
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logio(Tc) 



Figure 3: Relation between particle mass and central temperature. This graph displays 
the relation between logio(T'c) (in K) and logio(n^) (in keV’s) that follows from equation 
( b.4|) for a terminal velocity vq = 200 km/sec. Almost identical plots are obtained for 
other velocities in the observed range 125 km/sec < Vq < 250 km/sec. The circle and box 
symbols respectively denote the proton and electron mass yielding central temperatures 
of the order Tc ~ 10®, 10^ K. The central temperature for light particles in the range 
0.5keV < m < 100keV is less than 300 K (rectangle in the left), while for massive 
sypersymmetric particles in the range 1 GeV <m< 100 GeV, we have as large as 10® 
K (rectangle on the right). However, such high temperatures cannot rule out these weakly 
interactive particles as components of the dark matter MB gas. 
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